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The last decade has witnessed significant advances in our observational un- 
derstanding of the earliest stages of low-mass star formation. The advent of 
sensitive receivers on large radio telescopes such as the JCMT and IRAM 
30m MRT has led to the identification of young protostars at the beginning 
of the main accretion phase ('Class 0' objects), and has made it possible to 
probe, for the first time, the inner density structure of pre-collapse cores. 
Class objects are characterized by strong, centrally-condensed dust contin- 
uum emission at submillimeter wavelengths, very little emission shortward 
of ~ 10 fim, and powerful jet-like outflows. Direct evidence for gravitational 
infall has been observed toward several of them. They are interpreted as 
accreting protostars which have not yet accumulated the majority of their 
final stellar mass. In contrast to protostars, pre-stellar cores have flat inner 
density profiles, suggesting the initial conditions for fast protostellar collapse 
depart sometimes significantly from a singular isothermal sphere. In the case 
of non-singular initial conditions, the beginning of protostellar evolution is 
expected to feature a brief phase of vigorous accretion/ejection which may 
coincide with Class objects. In addition, submillimeter continuum imaging 
surveys of regions of multiple star formation such as Ophiuchus and Ser- 
pens suggest a picture according to which each star in an embedded cluster 
is built from a finite reservoir of mass and the associated IMF is primarily 
determined at the pre-stellar stage of evolution. 



I. INTRODUCTION 

The formation of low-mass stars is believed to involve a series of 
conceptually different stages (e.g. Larson 1969, Shu, Adams, & Lizano 
1987). The first stage corresponds to the fragmentation of a molecular 
cloud into a number of gravitationally-bound cores which are initially 
supported against gravity by a combination of thermal, magnetic, and 
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turbulent pressures (e.g. Mouschovias 1991, Shu et al. 1987). These 
pre-stellar condensations/fragments form and evolve as a result of a 
still poorly understood mechanism, involving ambipolar diffusion (e.g. 
Mouschovias 1991), the dissipation of turbulence (e.g. Nakano 1998), 
and/or an outside impulse (e.g. Bonnell ct al. 1997). Once such a 
condensation becomes gravitationally unstable and collapses, the main 
theoretical features of the ensuing dynamical evolution have been known 
since the pioneering work of Larson (1969). During a probably brief 
initial phase, the released gravitational energy is freely radiated away 
and the collapsing fragment stays roughly isothermal. This "runaway" 
isothermal collapse phase tends to produce a strong central concentra- 
tion of matter with a radial density gradient approaching p oc r~ 2 at 
small radii essentially independently of initial conditions (e.g. Whit- 
worth & Summers 1985, Blottiau et al. 1988, Foster & Chevalier 1993). 
It ends with the formation of an opaque, hydrostatic protostellar ob- 
ject in the center (e.g. Larson 1969, Boss & Yorke 1995, Bate 1998). 
Numerical simulations in fact predict the successive formations of two 
hydrostatic objects, before and after the dissociation of molecular hy- 
drogen respectively (see Boss & Yorke 1995), but we will not distinguish 
between them here. One then enters the main accretion phase during 
which the central object builds up its mass (M*) from a surrounding in- 
falling envelope (of mass M env ) and accretion disk, while progressively 
warming up. In this chapter, we will refer to the system consisting of 
the central object, plus envelope and disk as an accreting protostar. 
The youngest accreting protostars have M onv >> M*, and radiate the 
accretion luminosity L acc w GM*M acc /i?*. In the 'standard' theory 
of isolated star formation (Shu et al. 1987, 1993), the collapse ini- 
tial conditions are taken to be (static) singular isothermal spheroids 
(p ~ {a 2 /2irG)r~ 2 , cf. Li & Shu 1996, 1997), there is no runaway col- 
lapse phase, and the accretion rate M acc is constant in time ~ a 3 /G, 
where o is the effective isothermal sound speed. With other collapse 
initial conditions, the accretion rate is generally time-dependent (see 
III-D below). 

Observations have shown that the main accretion phase is always 
accompanied by a powerful ejection of a small fraction of the accreted 
material in the form of prominent bipolar jets/outflows (e.g. Bachillcr 
1996). These outflows are believed to carry away the excess angular 
momentum of the infalling matter (e.g. Konigl & Pudritz, this vol- 
ume). When the central object has accumulated most ( ^ 90%) of its 
final, main-sequence mass, it becomes a pre-main sequence (PMS) star, 
which evolves approximately at fixed mass on the Kelvin-Helmholtz 
contraction timescale (e.g. Stahler & Walter 1993). (Note that, during 
the protostellar accretion phase, stars more massive than a few 0.1 Mq 
start burning deuterium, while stars with masses in excess of ~ 8 Mq 
begin to burn hydrogen - see Palla & Stahler 1991.) 
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The details of the earliest stages outlined above are still poorly 
known. Improving our understanding of these early stages is of prime 
importance since to some extent they must govern the origin of the 
stellar initial mass function (IMF). 

Observationally, it is by comparing the structure of starless dense 
cores with that of the envelopes surrounding the youngest stellar ob- 
jects that one may hope to estimate the initial conditions for proto- 
stellar collapse. The purpose of this chapter is to review several major 
advances made in this field over the last decade, thanks mostly to 
ground-based (sub)millimcter continuum observations. We discuss re- 
sults obtained on pre-stellar cores and young accreting protostars in 
Sect. II and Sect. Ill, respectively. We then combine these two sets of 
results and conclude in Sect. IV. 

II. PRE-STELLAR CORES 

A. Definition and Identification 

The pre-stellar stage of star formation may be defined as the phase 
in which a gravitationally bound core has formed in a molecular cloud, 
and evolves toward higher degrees of central condensation, but no cen- 
tral hydrostatic protostellar object exists yet within the core. 

A pioneering survey of isolated dense cores in dark clouds was 
carried out in transitions of NH 3 by Myers and co-workers (see Myers 
& Benson 1983, Benson & Myers 1989 and references therein), who 
catalogued about 90 cores. These were separated into starless cores and 
cores with stars (Bcichman et al. 1986), on the basis of the presence 
or absence of an embedded source detected by IRAS. The starless NH3 
cores were identified by Beichman et al. as the potential sites of future 
isolated low-mass star formation. Other dense core surveys have been 
carried out by Clemens & Barvainis (1988), Wood et al (1994), Bourke 
et al (1995a,b), Lee & Myers (1999), and Jessop & Ward-Thompson 
(1999). 

Using the 15 m James Clerk Maxwell Telescope (JCMT), Ward- 
Thompson et al. (1994) observed the 800-^m dust continuum emission 
from about 20 starless NH 3 cores from the Benson & Myers list, map- 
ping 4 of the cores, and showed that they have larger FWHM sizes 
than, but comparable masses to the envelopes of the youngest proto- 
stars (Class sources - see III below). This is consistent with starless 
NH3 cores being pre-stellar in nature and the precursors of protostars 
(see also Mizuno et al. 1994). Ward-Thompson et al. also demon- 
strated that pre-stellar cores do not have density profiles which can 
be modelled by a single scale-free power law, but instead have fiat in- 
ner radial density profiles, suggestive of magnetically-supported cores 
contracting by ambipolar diffusion (see Mouschovias 1991, 1995 and 
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references therein). Recent molecular line spectroscopy of several pre- 
stellar cores (e.g. Tafalla et al. 1998 and Myers et al. this volume) 
appears to support the argument that they are contracting, but more 
slowly than the infall seen toward Class protostars (e.g. Mardones et 
al. 1997- see III-C). 

The 800 ,[mi study by Ward-Thompson et al. (1994) also suggests 
that wide-field submillimcter continuum imaging may be a powerful 
tool to search for new pre-stellar cores in the future (cf. Ristorcelli et 
al. 1998). 

B. Spectral Energy Distributions and Temperatures 

The advent of the Infrared Space Observatory (ISO) of the Euro- 
pean Space Agency (ESA), has allowed pre-stellar cores to be studied 
in the far-infrared for the first time (Ward-Thompson et al. 1999a 
in prep), since these cores were not detected by IRAS. Likewise, the 
Submillimetre Common User Bolometer Array (SCUBA) camera on 
the JCMT has allowed pre-stellar cores to be observed in the submm 
with greater signal to noise than ever before (Ward-Thompson et al. 
1999b in prep). Pre-stellar cores emit almost all of their radiation in the 
FIR/submm/mm regimes, so the combination of these two instruments 
provides a unique opportunity to study them. 

As an illustration, Figure 1 shows a series of images of the L1544 
pre-stellar core at 90 & 200 ^m (from ISO), 850 ^m (from SCUBA), and 
1.3 mm (from IRAM). The core is clearly detected at 200-1300 ^m, but 
is almost undetected at 90 /im. This shows that the core is very cold, 
and its dust temperature can be obtained from fitting a modified black- 
body to the observed emission. The spectral energy distribution (SED) 
of LI 544 in the far-infrared and submillimcter wavelength regimes is 
shown in Fig. 2. The solid line is a grey-body curve of the form: 

S v = B u {T dust ) [1 - exp(-r„)] Aft , 

where B v (Tdust) is the Planck function at frequency v for a dust tem- 
perature Tdust , t v — k v S is the dust optical depth through a mass 
column density £ (see below), and AO is the source solid angle. In this 
simple modelling, the dust opacity per unit (gas + dust) mass column 
density, k u , is assumed to scale as with f3 — 1.5-2 as usually appro- 
priate in the submillimeter range (e.g. Hildcbrand 1983). For L1544, a 
good fit to the SED is obtained with Tdust = 13 K and (3 = 2. Similar 
results are obtained in other starless cores. This confirms the lack of 
any warm dust in such cores, and consequently the lack of any embed- 
ded protostcllar object. The (sub)millimeter data show a morphology 
similar to the ISO images, at much higher resolution, indicating that 
the same dust is being observed at all wavelengths. Consequently the 
temperature derived from the SED is representative of all the emitting 
dust and can be used to convert submillimcter fluxes into estimates of 
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dust masses, and hence to gas masses. 

C. Mass and Density Structure 

Dust emission is generally optically thin at (sub) millimeter wave- 
lengths, and hence is a direct tracer of the mass content of molecular 
cloud cores. For an isothermal dust source, the total (gas + dust) mass 
M(r < R) contained within a radius R from the center is related to the 
submillimeter flux density S u (9) integrated over a circle of projected 
angular radius 6 = R/d by: 

M(r < R) = ttR 2 < E > R = [S v (0) d 2 ]/^ B v (T dust )], 

where < £ >_r is the average mass column density. The dust opacity n u 
is somewhat uncertain, but the uncertainties are much reduced when 
appropriate dust models are used (see Henning et al. 1995 for a review). 
For pre-stellar cores of intermediate densities (n# 2 ^ 10 5 cm -3 ), n v is 
believed to be close to K\.z = 0.005 cm 2 g _1 at 1.3 mm (e.g. Hildebrand 
1983, Preibisch et al. 1993). In denser cloud cores and protostellar 
envelopes, grain coagulation and the formation of ice mantles make k v 
a factor of <~ 2 larger, i.e., K1.3 = 0.01 cm 2 g _1 assuming a gas-to-dust 
mass ratio of 100 (e.g. Ossenkopf & Henning 1994). A still higher 
value, /C1.3 = 0.02 cm 2 g _1 , is recommended in protoplanetary disks 
(Beckwith et al. 1990, Pollack et al. 1994). 

Following this method, FWHM masses ranging from ~ 0.5 M 
to ~ 35 Mq are derived for the 9 isolated cores mapped by Ward- 
Thompson, Motte, & Andre (1999 - hereafter WMA99). 

Although observed pre-stellar cores are generally not circularly or 
elliptically symmetric (see, e.g., Fig. 1), one can still usefully constrain 
their radial density profiles by averaging the (sub)millimeter emission 
in circular or elliptical annuli. 

Figure 3a shows the azimuthally-averaged radial intensity profile 
of the pre-stellar core L1689B at 1.3 mm, compared to the profile of a 
spherical isothermal core model with p(r) oc r~ 2 . The model intensity 
profile results from a complete simulation taking into account both the 
observing technique (dual-beam mapping) and the reduction method 
(cf. Andre, Ward-Thompson, & Motte 1996). We see that L1689B 
exhibits the familiar radial profile of pre-stellar cores, with a flat in- 
ner region, steepening toward the edges (Ward-Thompson et al. 1994, 
Andre et al. 1996, WMA99). In this representative example, the radial 
density profile inferred assuming a constant dust temperature is as flat 
as p(r) cx r~ 0A (if the 3-D core shape is disk-like) or p(r) oc r -12 (if 
the core shape is spheroidal) at radii less than Rfi at ~ 4000 AU, and 
approaches p(r) oc r -2 only between ~ 4000 AU and ~ 15000 AU. 

More recently, it has been possible to constrain the outer density 
gradient of starless cores through absorption studies in the mid-infrared 
with ISOCAM (e.g. Abergeletal. 1996, 1998, Bacmann ct al. 1998). It 
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appears that isolated pre-stellar cores are often characterized by sharp 
edges, steeper than p oc r~ 3 or p oc r~ 4 , at radii R £ 15000 AU. 

These features of pre-stellar density structure, i.e., inner flatten- 
ing and sharp outer edge, are qualitatively consistent with models 
of magnetically-supported cores evolving through ambipolar diffusion 
prior to protostar formation (e.g. Ciolek & Mouschovias 1994, Basu & 
Mouschovias 1995), although the models generally require fairly strong 
magnetic fields (<~ 100 pG). Alternatively, the observed structure may 
also be explained by models of thermally supported self-gravitating 
cores interacting with an external UV radiation field (e.g. Falgarone & 
Puget 1985, Chieze & Pincau dcs Forets 1987). 

D. Lifetimes 

Beichman et al. (1986) used the ratio of numbers of starless cores 
to numbers of cores with embedded IRAS sources to estimate their 
relative timescales. They found roughly equal numbers of cores with 
and without IRAS sources. They estimated the lifetime of the stage of 
cores with stars, based on T Tauri lifetimes and prc-main sequence HR 
diagram tracks (e.g. Stahler 1988). Based on this, they estimated the 
lifetime of the pre-stellar core phase to be a few times 10 6 yr. 

However, the lifetime of a pre-stellar core depends on its central 
density. Figure 4 (taken from Jessop & Ward-Thompson 1999) shows 
the estimated lifetime of starless cores for each of six dark cloud surveys 
mentioned in II-A above, versus the mean volume density of cores in 
each sample. The lifespan of cores without stars in each sample was 
estimated from the fraction of cores with IRAS sources, using the same 
method as Beichman ct al. (1986). An anti-correlation between lifetime 
and density is clearly apparent in Fig. 4. The solid line has the form 
t oc p~°' 75 , while the dashed line is of the form t oc p~ - 5 . These 
two forms are expected for cores evolving on the ambipolar diffusion 
timcscale tAD °^ x e (where x e is the ionisation fraction - e.g. Nakano 
1984), if the dominant ionisation mechanism is cosmic ray ionisation or 
UV ionisation, respectively (e.g. McKee 1989). 

More details about the evolution of pre-stellar cores at higher den- 
sities can be inferred from the results of the submm continuum SCUBA 
survey of Ward-Thompson et al. (1999b). In this survey, 17 of the 38 
NH3 cores without IRAS sources from Benson & Myers (1989) were 
detected by SCUBA at 850^m. Ward-Thompson ct al. (1999b) esti- 
mate that the 17 SCUBA detections all have central densities between 
~10 5 cm -3 and ~10 6 cm~ 3 , whilst the 21 non-detections must have 
lower central densities, typically between ~10 4 cm~ 3 and ^10 5 cm~ 3 
(see Benson & Myers 1989, Butner et al. 1995, and references therein). 
Consequently, they deduce that the lifetime of these two phases - cen- 
tral density increasing from ~10 cm~ 3 to ~10 cm compared to cen- 
tral density of ^-40 5 cm~ 3 until the formation of a protostellar object 
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at the center - must be roughly equal. This can be compared with the 
predictions of ambipolar diffusion models. 

Figure 3b is taken from Ciolek & Mouschovias (1994) and shows 
the radial density profile predicted by an ambipolar diffusion model at 
different evolutionary stages (to to t 6 ). The stage at which the central 
density is ~10 4 cm -3 corresponds to time t\ and the stage at which 
the central density is ~10 5 cm -3 corresponds to time t2- The time at 
which a protostellar object forms is effectively t&. In this model the 
time taken to go from t\ to ti (~2x 10 6 yr) is six times longer than 
the time taken to go from ti to t§. Some discrepancy could perhaps 
be accounted for by the statistical errors associated with our source 
number counting technique, but the ratio between the two timescales 
should be fairly robust: The model predicts that SCUBA should only 
have detected ~l/7 of the cores, whereas it detected half of the sample. 

We are left with the conclusion that cores at central densities of 
order ~ 1(P cm -3, evolve more slowly than ambipolar diffusion models 
predict - i.e. the cores experience more support than a simple static 
magnetic field can provide. The extra support could perhaps be pro- 
vided by turbulence which generates non-static magnetic fields (e.g. 
Gammic & Ostriker 1996, Nakano 1998, Balsara et al. 1998). 

E. Pre-Stellar Condensations in Star-Forming Clusters 

In regions of multiple star formation, submillimctcr dust contin- 
uum mapping has revealed a wealth of small-scale cloud fragments, 
sometimes organized along filaments (e.g. Mezger et al. 1992, AWB93, 
Casali et al. 1993, Launhardt et al. 1996, Chini et al. 1997b, John- 
stone & Bally 1999). Such fragmentation along filaments has not been 
observed in Taurus, but examples do exist in young embedded clusters 
forming primarily low-mass stars like p Ophiuchi (Motte, Andre, & 
Neri 1998). The individual fragments, which are denser (< n > £ 10 6 - 
10 7 cm -3 ) and more compact (a few 1000 AU in size) than the isolated 
pre-stellar cores discussed above, often remained totally undetected (in 
emission) by IRAS or ISO in the mid- to far-IR. Since molecules tend 
to freeze out onto dust grains at low temperatures and high densities, 
(sub)millimeter dust emission may be the most effective tracer of such 
condensations (e.g. Mauersberger et al. 1992). 

The most centrally-condensed of these starless fragments have been 
claimed to be isothermal protostars, i.e., collapsing condensations with 
no central hydrostatic object (see I) (e.g. Mezger et al. 1992, Launhardt 
et al. 1996, Motte et al. 1998). Good examples are FIR 3 and FIR 4 
in NGC 2024, OphA-SMl and OphE-MM3 in p Oph, LBS 17-SM in 
Orion B, or MMS1 and MMS4 in OMC-3. This isothermal protostar 
interpretation remains to be confirmed, however, by observations of 
appropriate spectral line signatures (cf. Myers et al. this volume). 

Furthermore, the advent of large-format bolometer arrays now 



8 



Andre, Ward-Thompson, Barsony 



makes possible systematic studies of the genetic link between prc-stcllar 
cloud fragments and young stars. Figure 5 is a 1.3 mm continuum wide- 
field mosaic of the p Oph cloud (Motte et al. 1998) showing a total 
of 100 structures with characteristic angular scales of ~ 15"-30" (i.e., 
~ 2000-4000 AU) , which are associated with 59 starless condensations 
(undetected by ISO in the mid-IR) and 41 embedded YSOs (detected 
at IR or radio continuum wavelengths). 

Comparison of the masses derived from the 1.3 mm continuum 
(from ^0.05M Q to ^3M Q ) with Jeans masses suggests that most of the 
59 starless fragments are close to gravitational virial equilibrium with 
M/M vir £ 0.3-0.5 and will form stars in the near future. These pre- 
stcllar condensations generally have flat inner density profiles like iso- 
lated pre-stellar cores, but are distinguished by compact, finite sizes of 
a few thousand AU. The typical fragmentation lengthscale derived from 
the average projected separation between condensations is ~ 6000AU 
in p Oph. This is ~ 5 times smaller than the radial extent of isolated 
dense cores in the Taurus cloud (see Gomez et al. 1993). 

Figure 6 shows the mass distribution of the 59 p Oph pre-stellar 
fragments. It follows approximately AN / AM cx M -1,5 below ~0.5M Q , 
which is similar to the clump mass spectrum found by large-scale molec- 
ular line studies (e.g. Blitz 1993 and Williams et al. this volume). 
The novel feature, however, is that the fragment mass spectrum found 
at 1.3 mm in p Oph appears to steepen to AN/ AM cx M~ 2 - 5 above 
^0.5M Q . A similarly steep mass spectrum above ~O.5M was obtained 
by Testi & Sargent (1998) for compact 3 mm starless condensations in 
the Serpens core. These pre-stellar mass spectra resemble the shape 
of the stellar initial mass function (IMF), which is known to approach 
AAV AM* « M ~ 2 - 7 for 1 M Q < M* < 10 M Q and AN/ AM* cx M" 1 - 2 
for 0.1 M Q < M* < 1 M Q (e.g. Kroupa et al. 1993, Tinncy 1993, 1995, 
Scalo 1998, see also Meyer et al. this volume). Given the factor of ~ 2 
uncertainty on the measured pre-stellar masses, such a resemblance is 
remarkable and suggests that the IMF of embedded clusters is primar- 
ily determined at the pre-stellar stage of star formation (see also IV 
below). 

III. THE YOUNGEST PROTOSTARS 

A. Class Protostars and Other YSO Stages 

1. Infrared YSO classes. In the near-/mid-infrared, three broad 
classes of young stellar objects (YSOs) can be distinguished based on 
the slope am = dlog(AF A )/dlog(A) of their SEDs between 2.2 pm and 
10-25 /im, which are interpreted in terms of an evolutionary sequence 
(Lada & Wilking 1984, Lada 1987). Going backward in time, Class III 
(aiR < —1.5) and Class II (—1.5 < aiR < 0) sources correspond to 
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PMS stars ("Weak" and "Classical" T Tauri stars, respectively) sur- 
rounded by a circumstellar disk (optically thin and optically thick at 
A ^ 10 yum, respectively), but lacking a dense circumstellar envelope 
(see Andre & Montmerle 1994 - hereafter AM94). The youngest YSOs 
detected at 2 /xm are the Class I sources, which are characterized by 
Qir > (e.g. Wilking, Lada, & Young 1989 - WLY89), and the close 
association with dense molecular gas (e.g. Myers et al. 1987). Class I 
objects are now interpreted as relatively evolved protostars with typical 
ages ~ 1-2 x 10 5 yr (e.g. Barsony & Kenyon 1992, Greene et al. 1994, 
Kcnyon & Hartmann 1995), surrounded by both a disk and a diffuse 
circumstellar envelope of substellar ( £ 0.1-0.3 M Q ) mass (Whitney 
& Hartmann 1993, Kenyon et al. 1993b, AM94, Lucas & Roche 1997). 
Their SEDs are successfully modeled in the framework of the "stan- 
dard" theory of isolated protostars (e.g. Adams, Lada, & Shu 1987, 
Kenyon et al. 1993a), in agreement with the idea that they derive a 
substantial fraction of their luminosity from accretion (see also Greene 
& Lada 1996 and Kenyon et al. 1998). 

2. Class protostars. Several condensations detected in submil- 
limeter dust continuum maps of molecular clouds (e.g. II-E) appear 
to be associated with formed, hydrostatic YSOs and have been desig- 
nated "Class 0" protostars (Andre, Ward-Thompson, & Barsony 1993 
- AWB93). Specifically, Class objects are defined by the following 
observational properties (AWB93): 

• (i) Indirect evidence for a central YSO, as indicated by, e.g., the 
detection of a compact centimeter radio continuum source, a collimated 
CO outflow, or an internal heating source. 

• (ii) Centrally peaked but extended submillimeter continuum 
emission tracing the presence of a spheroidal circumstellar dust en- 
velope (as opposed to just a disk). 

• (iii) High ratio of submillimeter to bolometric luminosity suggest- 
ing the envelope mass exceeds the central stellar mass: \j S mml^boi > 
0.5%, where L 5mm is measured longward of 350 /im. In practice, this of- 
ten means a SED resembling a single temperature blackbody at T ~ 15- 
30 K (see Fig. 2). 

Property (i) distinguishes Class objects from the pre-stellar cores 
and condensations discussed in Sect. II. In particular, deep VLA ob- 
servations reveal no compact radio continuum sources in the centers of 
pre-stellar cores (Bontemps 1996, Yun et al. 1996). Properties (ii) and 
(iii) distinguish Class objects from more evolved (Class I and Class II) 
YSOs. As shown by AWB93, the \j smm /^boi ratio should roughly track 
the ratio M env /M+ of envelope to stellar mass, and may be used as 
an evolutionary indicator (decreasing with time), for low-luminosity 
(Lboi ^, 50 Lq) embedded YSOs. Criterion (iii) approximately selects 
objects which have M env /M* > 1, assuming plausible relations be- 
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tween and M* on the one hand and between L smm and M env on 
the other hand (see AWB93 and AM94). [A roughly equivalent crite- 
rion is M env /L bo i > 0.1 Mq/Lq.] Class objects are therefore excellent 
candidates for being very young accreting protostars in which a hydro- 
static core has formed but not yet accumulated the majority of its final 
mass. In practice, most of the confirmed Class objects listed in Ta- 
ble 1 have L smm /Lh ; >> 0.5% and are likely to be at the beginning of 
the main accretion phase with M env » M* (see Fig. 7b). 

3. Evolutionary diagrams for embedded YSOs. Combining infrared 
and submillimeter data, it is therefore possible to define a complete, em- 
pirical evolutionary sequence (Class — > Class I — > Class II — > Class III) 
for low-mass YSOs, which likely correspond to conceptually different 
stages of evolution: (early) main accretion phase, late accretion phase, 
PMS stars with protoplanetary disks, PMS stars with debris disks (see 
AM94). This sequence is quasi-continuous and may be parameterized 
by the "bolometric temperature" , T^oi, defined by Myers & Ladd (1993) 
as the temperature of a blackbody having the same mean frequency as 
the observed YSO SED. Myers & Ladd proposed to use the Lboi-Tboi 
diagram for embedded YSOs as a direct analog to the H-R diagram 
for optically visible stars. As shown by Chen et al. (1995, 1997), YSOs 
with known classes have distinct ranges of Tf, ;: < 70 K for Class 0, 70- 
650 K for Class I, 650-2880 K for Class II, and > 2880 K for Class III 
(e.g. Fig. 7a). The evolution of Tt, i and L\, i from the Class stage to 
the zero-age main sequence has been modelled in the context of various 
envelope-dissipation scenarios by Myers et al. (1998). 

A perhaps more direct approach to tracking the circumstellar evo- 
lution of YSOs is to use the circumstellar mass M c * derived from 
(sub) millimeter continuum measurements of optically thin dust emis- 
sion. Such measurements show that M c * (= M env + Mdisk) is generally 
dominated by M env in Class 0/Class I sources (e.g. Terebey et al. 1993) 
and decreases by a factor ~ 5-10 on average from one YSO class to the 
next (AM94). In the spirit of the L smm /Lb i evolutionary indicator of 
AWB93, Saraceno et al. (1996a) proposed the L smm -Lb i (or equiva- 
lcntly M env ~Lboi) diagram as an alternative evolutionary diagram for 
self-embedded YSOs. While L smm and M env are well correlated with 
Lboi for the majority of embedded YSOs (e.g. Reipurth et al. 1993), 
Class objects clearly stand out from Class I sources in this diagram as 
objects with excess (sub)millimeter emission, i.e., excess circumstellar 
material (see Fig. 7b). Moreover, one may compare the locations of 
observed embedded sources in the M env -Lboi diagram with simple pro- 
tostellar evolutionary tracks (Saraceno et al. 1996a, b). Qualitatively 
at least, scenarios in which the mass-accretion rate decreases with time 
for a given protostar (Bontcmps et al. 1996a, Myers et al. 1998 - see 
also III-D below) and increases with the mass of the initial pre-collapse 
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fragment (e.g. Myers & Fuller 1993, Rcipurth ct al. 1993, Saraceno 
et al. 1996b) yield tracks in better agreement with observations than 
the constant-rate scenario discussed by Saraceno et al. (1996a). In 
particular, the peak accretion luminosity is reduced by a factor ~ 2-4 
compared to the constant-rate scenario (cf. Bontemps et al. 1996a, 
Myers et al. 1998, and Fig. 7b), which agrees better with observed 
luminosities (e.g. Kenyon & Hartmann 1995). 

Inclination effects may a priori affect the positions of individual 
protostellar objects in these evolutionary diagrams. In particular, it has 
been claimed that some Class I sources may potentially look like Class 
objects when observed at high inclination angles to the line of sight (e.g. 
Yorke, Bodenhcimcr, & Laughlin 1995, Sonnhalter, Preibisch, & Yorkc 
1995, Men'shchikov & Henning 1997). However, the fact that Class 
objects are associated with an order of magnitude more powerful out- 
flows than Class I sources (see III-D below) confirm that these two 
types of YSOs differ qualitatively from each other. Furthermore, some 
Class sources are known to have small inclination angles (e.g. Cabrit 
& Bcrtout 1992, Greaves et al. 1997, Wolf-Chase et al. 1998). We 
also stress that existing (sub)millimeter maps of dust continuum and 
C 18 emission provide direct evidence that both Class I and Class 
objects are self-embedded in substantial amounts of circumstcllar mate- 
rial distributed in spatially resolved, spheroidal envelopes (e.g. AM94, 
Chen et al. 1997, Ladd et al. 1998, Dent et al. 1998 - see also III-B 
below). This material has the ability to absorb the optical and near- 
IR emission from the underlying star/disk system and to reradiate it 
quasi-isotropically at longer far-IR and (sub)mm wavelengths. In such 
a self-embedded configuration, viewing-angle effects are minimized, as 
confirmed by radiative transfer calculations (e.g. Efstathiou & Rowan- 
Robinson 1991, Yorke et al. 1995). Physically, this is because the 
bulk of the luminosity emerges at long wavelengths where most of the 
emission is effectively optically thin. However, the short-wavelength 
emission from the inner star/disk remains very dependent on viewing 
angle, implying that T^oi estimates should be somewhat more sensitive 
to orientation effects than h^oi and L smm /Lfc ;. 

4- Protostar surveys and lifetime estimates. Based on the key 
attributes of Class protostars (see 2. above), various strategies can 
be used to search for and discover new candidates: (sub)millimeter 
continuum mapping (e.g. Mezgeretal. 1992, Casali et al. 1993, Sandcll 
& Weintraub 1994, Reipurth et al. 1996, Launhardt ct al. 1996, Chini 
et al. 1997ab, Motte et al. 1998, Andre et al. 1999), HIRES processing 
of the IRAS data (Hurt & Barsony 1996, O'Linger et al. 1999), deep 
radio continuum VLA surveys (e.g. Leous et al. 1991, Bontemps et al. 
1995, Yun ct al. 1996, Moreira et al. 1997, Gibb 1999), CO mapping 
(e.g. Bachiller et al. 1990, Andre et al. 1990, Bourke et al. 1997), and 
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large-scale near-IR/optical imaging of shocked H2 and [SII] emission 
(e.g. Hodapp & Ladd 1995, Yu, Bally, & Devine 1997, Wilking et al. 
1997, Gomez et al. 1998, Stanke et al. 1998, Phelps & Barsony 1999). 

Class objects appear to be short-lived compared to both pre- 
stellar cores/fragments (see Sect. II above) and Class I near-IR sources. 
For example, in the p Oph main cloud where the IRAM 30 m mapping 
study of Motte et al. (1998) provides a reasonably complete census 
of both pre- and proto-stellar condensations down to ^0.1 M Q , there 
are only two good Class candidates (including the prototypical object 
VLA 1623 - AWB93), while there are ~ 15-30 near-IR Class I sources 
(e.g. WLY89, AM94, Greene et al. 1994, Barsony & Ressler 1999, 
Bontemps et al. 1999). 

Under the assumption that p Oph is representative and forming 
stars at a roughly constant rate, the lifetime of Class objects should 
thus be approximately an order of magnitude shorter than the lifetime 
of Class I sources (see 1. above), i.e., typically ~ 1-3 xlO 4 yr. The 
jet-like morphology and short dynamical timescales of Class outflows 
are consistent with this estimate (e.g. Barsony et al. 1998, see III-D 
below). A lifetime as short as a few 10 4 yr supports the interpretation 
of Class objects as very young accreting protostars (see, e.g., Fletcher 
& Stahler 1994 and Barsony 1994). 

Regions like the Serpens, Orion, or Perseus/NGC1333 complexes 
seem to be particularly rich in Class objects, and this is probably 
indicative of fairly high levels of ongoing, likely induced star formation 
activity (e.g. Hurt & Barsony 1996, Chini et al. 1997b, Yu et al. 
1997, Barsony et al. 1998). For instance, we estimate the current star 
formation rate in Orion-OMC-3 to be ~ 2 x 10~ 3 M Q yr -1 , which is 1-2 
orders of magnitude larger than the star formation rates characterizing 
the Trapezium, NGC 1333, and IC 348 near-IR clusters when averaged 
over £ 10 6 yr periods (see Lada, Alves, & Lada 1996). 

B. Density Structure of the Protostellar Environment 

1. Envelope. In contrast to pre-stellar cores, the envelopes of low- 
mass Class and Class I protostars are always found to be strongly 
centrally-condensed and do not exhibit any inner flattening in their 
(sub)millimeter continuum radial intensity profiles. In practice, this 
means that, when protostars are mapped with the resolution of the 
largest single-dish telescopes, the measured peak flux density is typi- 
cally a fraction £ 20% of the flux integrated over five beam widths. 
For comparison, the same fraction is H 10% for pre-stellar cores. 
(Sub)millimeter continuum maps indicate that protostellar envelopes 
in regions of isolated star formation such as Taurus have radial density 
gradients generally consistent with p(r) oc r~ p withp ~ 1.5-2 over more 
than ~ 10000-15000 AU in radius (e.g. Walker, Adams, & Lada 1990, 
Ladd et al. 1991, Motte 1998, Motte et al. 1999, see also Mundy et 
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al. this volume) . The estimated density gradient thus agrees with most 
collapse models which predict a value of p between 1.5 and 2 during the 
protostellar accretion phase (e.g. Whitworth & Summers 1985). Some 
studies have, however, inferred shallower density gradients (p ~ 0.5-1 - 
e.g. Barsony & Chandler 1993, Chandler, Barsony, & Moore 1998). In 
any case, the densities and masses measured for the envelopes around 
the bona-fide Class I objects of Taurus appear to be consistent within 
a factor of ~ 4 with the predictions of the "standard" inside-out col- 
lapse theory (e.g. Shu et al. 1993) for ~ 10 5 yr-old, isolated protostars 
(Motte 1998, Motte et al. 1999). 

The situation is markedly different in star-forming clusters. In 
p Ophiuchi in particular, the circumstellar envelopes of Class I and 
Class protostars are observed to be very compact: they merge with 
dense cores, other envelopes, and/or the diffuse ambient cloud at a. finite 
radius R out < 5 000 AU (Motte et al. 1998). This is > 3 times smaller 
than the collapse expansion wavefront at a 'Class I age' of ~ 2 x 10 5 yr 
in the standard theory of isolated protostars, emphasizing the fact that 
each YSO has a finite 'sphere of influence' in p Oph. Similar results 
were obtained in the case of the Perseus Class sources NGC1333-4A, 
NGC1333-2, L1448-C, and L1448-N (e.g. Motte 1998). Moreover, the 
envelopes of these Perseus protostars are 3 to 10 times denser than the 
singular isothermal sphere for a sound speed a = 0.2 kms -1 . This 
suggests that, prior to collapse, the main support against gravity was 
turbulent and/or magnetic in origin rather than purely thermal (see 
also Mardones et al. 1997). 

2. Disks and multiplicity. Many Class protostars are in fact 
multiple systems, when viewed at sub-arcsecond resolution, sharing a 
common envelope and sometimes a circumbinary disk (e.g. Looney et 
al. 1999 - see also col. 11 of Table 1 and chapter by Mundy, Looney, 
& Welch). These protobinaries probably formed by dynamical frag- 
mentation during (or at the end of) the isothermal collapse phase (e.g. 
Chapman et al. 1992, Bonnell 1994, Boss & Myhill 1995). Interest- 
ingly enough, only cores with inner density profiles as flat as p oc r _1 
or flatter (like observed pre-stellar cores - see Sect. II), can apparently 
fragment during collapse (Myhill & Kaula 1992, Boss 1995, Burkert et 
al. 1997). 

Despite the difficulty of discriminating between the disk and enve- 
lope components, existing (sub)millimeter continuum interferometric 
measurements suggest that the "disks" of Class objects are a factor 
of £ 10 less massive than their surrounding circumstellar envelopes 
(e.g. Chandler et al. 1995, Pudritz et al. 1996, Looney et al. 1999, 
Hogerheijde 1998, Motte 1998, Wilner & Lay this volume). 



14 



Andre, Ward-Thompson, Barsony 



C. Direct Evidence for Infall 

Rather convincing spectroscopic signatures of gravitational infall have 
been reported toward several Class objects, confirming their proto- 
stellar nature (e.g. Walker et al. 1986, Zhou et al. 1993, Gregersen et 
al. 1997, Mardones et al. 1997, see also col. 10 of Table 1). Inward 
motions can be traced by optically thick molecular lines which should 
(locally) exhibit asymmetric self-absorbed profiles skewed to the blue 
(see Myers, Evans, & Ohashi, this volume). The interpretation is often 
complicated by the simultaneous presence of rotation and/or outflow 
(e.g. Menten et al. 1987, Walker et al. 1994, Cabrit et al. 1996). 

A comprehensive survey of a sample of 47 embedded YSOs in 
H 2 CO(2i2 — In) and CS(2-1) suggests that infall is more prominent in 
Class than in Class I sources (Mardones et al. 1997): In these transi- 
tions, infall asymmetries are detected toward 40-50 % of Class objects 
but less than 10 % of Class I sources. This is qualitatively consistent 
with a decline of infall/accretion rate with evolutionary stage (see III— 
D below). However, a more recent survey by Gregersen et al. (1999) 
using HCO + (3-2) finds no difference in the fraction of sources with 
"blue profiles" between Class and Class I sources. The HCO+(3-2) 
line is more optically thick than H 2 CO(2i 2 — In) and CS(2-1), hence a 
better tracer of infall at advanced stages. This result shows that some 
infall is still present at the Class I stage but remains consistent with a 
decline of the net accretion rate with time. The outflow is so broad in 
Class I sources that there often appears to be little transfer of mass to 
the inner ~ 2000 AU radius region around these objects (e.g. Fuller et 
al. 1995b, Cabrit et al. 1996, Brown & Chandler 1999). 

D. Decline of Outflow and Inflow with Time 

1. Evolution from Class to Class I. Most, if not all, Class 
protostars drive powerful, "jet-like" CO molecular outflows (see, e.g., 
Bachiller 1996 and Richer et al., this volume, for reviews). The me- 
chanical luminosities of these outflows are often of the same order as the 
bolometric luminosities of the central sources (e.g. Curiel et al. 1990, 
AWB93, Barsony et al. 1998). In contrast, while there is good evidence 
that some outflow activity exists throughout the accretion phase (e.g. 
Terebey et al. 1989, Parker et al. 1991, Bontcmps et al. 1996a), the 
CO outflows from Class I sources tend to be much less powerful and 
less collimated than those from Class objects. 

In an effort to quantify the evolution of mass loss during the pro- 
tostellar phase, Bontcmps et al. (1996a - hereafter BATC) obtained 
and analyzed a homogeneous set of CO (2-1) outflow data around a 
large sample of low-luminosity (L^oi < 50 L ), nearby (d < 450 pc) 
self-embedded YSOs. Their results show that Class objects lie an 
order of magnitude above the well-known (e.g. Cabrit & Bertout 1992) 
correlation between outflow momentum flux (Fco) and bolometric lu- 
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minosity (Lbol) holding for Class I sources (see Fqo ~Lboi diagram 
shown in Fig. 5 of BATC). Furthermore, BATC found that F C o was 
well correlated with M env in their entire sample (including both Class I 
and Class sources). The same correlation was noted independently on 
other source samples by Moriarty-Schieven et al. (1994), Hogerheijde 
et al. (1998), and Hcnning & Launhardt (1998). As argued by BATC, 
this new correlation is independent of the Fqo ^^boi correlation and 
most likely results from a progressive decrease of outflow power with 
time during the accretion phase. This is illustrated in the normalized 
Fco c/Lboi versus M env /L° 6 j diagram of Fig. 8, which should be essen- 
tially free of any luminosity effect. 

Since magneto-centrifugal accretion/ejection models of bipolar out- 
flows (e.g. Shu et al. 1994, Ferreira & Pclletier 1995, Fiege & Henrikscn 
1996, Ouyed & Pudritz 1997) predict a direct proportionality between 
accretion and ejection, BATC proposed that the decline of outflow 
power with evolutionary stage seen in Fig. 8 reflects a corresponding de- 
crease in the mass-accretion/infall rate. The results of BATC indicate 
that Mjet declines from <~ 1CP 6 M Q yr _1 for the youngest Class pro- 
tostars to ~ 2 x 10~ 8 M Q yr _1 for the most evolved Class I sources, sug- 
gesting a decrease in M acc from <~ 10~ 5 M Q yr _1 to ~ 2 x 1CP 7 M Q yr" 1 
if realistic jet model parameters are adopted {M- ]ct / M acc ~ 0.1-0.3, 
Vjot ~ 100 km s -1 ). These indirect estimates of M acc for Class and 
Class I protostars should only be taken as indicative of the true evo- 
lutionary trend. Nevertheless, it is interesting to note that they agree 
well with independent estimates of the rates of envelope dissipation 
based on circumstellar mass versus age arguments (Ward-Thompson 

1996, Ladd, Fuller, & Dcane 1998). 

As illustrated by the evolutionary tracks of Fig. 7b, a decline of 
-Macc with time does not imply a higher accretion luminosity for Class 
objects compared to Class I sources because the central stellar mass is 
smaller at the Class stage and the stellar radius is likely to be larger 
(see Henriksen, Andre, & Bontemps 1997 - hereafter HAB97). 

2. Link with the collapse initial conditions. The apparent decay of 
M acc from the Class to the Class I stage may be linked with the den- 
sity structure observed for pre-stellar cores/condensations (Sect. II). 
(Magncto)hydrodynamic collapse models predict a time-dependent ac- 
cretion history when the radial density profile at the onset of fast pro- 
tostellar collapse differs from p oc r~ 2 (e.g. Foster & Chevalier 1993 
- FC93, Tomisaka 1996, McLaughlin & Pudritz 1997, HAB97, Basu 

1997, Saner, McKee, & Stahler 1997, Li 1998, Ciolek & Konigl 1998). 
In particular, when starting from Bonnor-Ebert-like initial conditions 
resembling observed pre-stellar cores, these studies find that supersonic 
infall velocities develop prior to the formation of the central hydrostatic 
protostellar object at t = (see, e.g., FC93). Observationally, this 
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early collapse phase should correspond to 'isothermal protostars' (see 
II-E for possible examples). During the protostellar accretion phase 
(t > 0), because of the significant infall velocities achieved at t < 0, 
M acc is initially higher than the standard ~ a 3 /G value obtained for 
the inside-out collapse of a static singular isothermal sphere (Shu 1977, 
see also Sect. I). The accretion rate then converges toward the stan- 
dard value of the Shu solution, and declines again below a 3 /G at late 
times if the reservoir of mass is finite. By comparison with the rough 
estimates of M aC c given above, it is tempting to identify the short pe- 
riod of energetic accretion (M acc ~ 10 x a 3 /G) predicted by the mod- 
els just after point-mass formation with the observationally-defincd 
Class stage (HAB97, see Fig. 8). In this view, the more evolved 
Class I objects would correspond to the longer period of moderate ac- 
cretion/ejection when the accretion rate approaches the standard value 
(M acc < a 3 /G - 2 x 1(T 6 M Q yr" 1 for a cloud temperature of - 10 K). 
Using simple pressure-free analytical calculations (justified since the 
inflow becomes supersonic early on), HAB97 could indeed find a good 
overall fit to the empirical accretion history inferred by BATC on the 
basis of CO outflow observations (see Fig. 8). 

In the absence of magnetic fields, FC93 have shown that the 
timescale for convergence to the standard accretion rate of Shu (1977) 
depends on the radius of the flat inner core (Rfi a t) relative to the 
outer radius of the initial pre-collapse condensation (R ou t)- They found 
that a phase of constant ~ a 3 /G accretion rate is achieved only when 
Rout/ Rfiat ~ 20, typically after ~ 10 free-fall times of the flat inner 
region, for a period lasting ~ 15 free-fall times when R ou t/ Rfiat = 20 
and progressively longer as R ou t/ Rfiat increases. Since observations 
suggest Rout/ Rfiat >$• 3 in Ophiuchus and R ou t/ Rfiat ~ 15 in Taurus 
(see II-C, III-B, and Motte et al. 1998), one may expect a marked 
time-dependence of M acc in Ophiuchus but a reasonable agreement 
with the constant accretion rate of the self-similar theory of Shu et al. 
(1987, 1993) in Taurus. Indeed, HAB97 note that there is a much bet- 
ter continuity between Class and Class I protostars in Taurus than 
in Ophiuchus (see also Andre 1997). 

Finally, we stress that the absolute values of M acc in the Class and 
Class I stages are presently quite uncertain. An alternative interpreta- 
tion of the evolution seen in Fig. 8 is that Class protostars accrete at 
a rate roughly consistent with ~ a 3 /G, while most Class I sources are 
in a terminal accretion phase with M acc £ 0.1 x a 3 /G, resulting from 
the finite effective reservoir of mass available to each object in clusters 
(e.g. Motte et al. 1998 and III-B) and/or from the effects of outflows 
dispersing the envelope (e.g. Myers et al. 1998, Ladd et al. 1998). 
Distinguishing between this possibility and that advocated above will 
require direct measurements of the mass accretion rates. 
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IV. CONCLUSIONS AND IMPLICATIONS 

The observational studies discussed in Sect. II demonstrate that 
pre-stellar cores/fragments are characterized by flat inner radial den- 
sity gradients. This, in turn, suggests the initial conditions for fast pro- 
tostellar collapse are non-singular, i.e., the density profile at the onset 
of collapse is not infinitely centrally condensed. (Sub)millimcter obser- 
vations also set strong constraints on protostellar evolution (Sect. III). 
The fact that young (Class 0) protostars drive more powerful outflows 
than evolved (Class I) protostars suggests that the mass accretion rate 
M acc decreases by typically a factor of ~ 5-10 from the Class to the 
Class I stage (III-D). Such a decline in M acc during the protostellar 
accretion phase may be the direct result of a flattened initial density 
profile (see III-D). Based on these observational constraints, we sug- 
gest that most protostars form in a dynamical rather than quasi-static 
fashion. 

The results summarized in this chapter also have broader impli- 
cations concerning, e.g., the origin of the IMF. As pointed out in 
Sect. II-E, the pre-stellar condensations observed in regions of mul- 
tiple star formation such as p Ophiuchi arc finite-size structures, typ- 
ically a few 1000 AU in radius, which are clearly not scale- free. This 
favors a picture of star formation in clusters in which individual pro- 
tostellar collapse is initiated in compact dense clumps resulting from 
fragmentation and resembling more finite-size Bonnor-Ebert cloudlets 
than singular isothermal spheres. Such condensations may correspond 
to dense, low-ionization pockets decoupling themselves from the parent 
molecular cloud as a result of ambipolar diffusion and/or the dissipa- 
tion of turbulence (e.g. Mouschovias 1991, Myers 1998, Nakano 1998). 
They would thus be free to undergo Jeans-like gravitational instabili- 
ties and collapse under the influence of external disturbances, of which 
there are many types in regions of multiple star formation (e.g. Pringlc 
1989, Whitworth ct al. 1996, Motte et al. 1998, Barsony et al. 1998). 
By contrast, the low-density regions of the ambient cloud, being more 
ionized, would remain supported against collapse by static and/or tur- 
bulent magnetic fields. The typical separation between individual con- 
densations should be of order the Jeans length in the parent cloud/core, 
in rough agreement with observations (e.g. Motte et al. 1998). 

In this obscrvationally-driven scenario of fragmentation and col- 
lapse, stars are built from bounded fragments which represent finite 
reservoirs of mass. The star formation efficiency within these frag- 
ments is high: most of their 'initial' masses at the onset of collapse 
end up in stars. If this is true, it implies that the physical mecha- 
nisms responsible for the formation of pre-stellar cores/condensations 
in molecular clouds, such as turbulent fragmentation (e.g. Padoan, 
Nordlund, & Jones 1997), play a key role in determining the IMF of 
embedded clusters. Such a picture, which we favor for regions of mul- 
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tiplc star formation, is consistent with some theoretical scenarios of 
protocluster formation (e.g. Larson 1985, Klessen, Burkert, & Bate 
1998). It need not be universal, however, and is in fact unlikely to 
apply to regions of isolated star formation like the Taurus cloud. In 
these regions, protostars may accrete from larger reservoirs of mass, 
and feed-back processes such as stellar winds may be more important 
in limiting accretion and defining stellar masses (e.g. Shu et al. 1987, 
Silk 1995, Adams & Fatuzzo 1996, Velusamy & Langer 1998). 

With the advent of major new facilities at far-IR and submillimcter 
wavelengths, the next decade promises to be at least as rich in observa- 
tional discoveries as the past ten years. By combining the capabilities of 
space telescopes such as FIRST with those of large ground-based arrays 
such as the LSA/MMA, it will be possible to study the detailed physics 
of complete samples of young protostars and pre-collapse fragments, in 
a variety of star-forming clouds, and down to the brown-dwarf mass 
regime. This should tremendously improve our global understanding 
of the initial stages of star formation in the Galaxy. 
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FIGURE CAPTIONS 

Figure 1. Dust continuum images of L1544 at 90 /im (a) and 200 pm (b) from 
ISOPHOT, at 850 pm (c) from SCUBA, and at 1.3 mm (d) from IRAM 
30 m. A polarization E-vector, perpendicular to the B field, is overlaid 
on the 850 /jm image (c), as measured with the SCUBA polarimeter. The 
observed morphology is consistent with a magnetically-supported core that 
has flattened along the direction of the mean magnetic field. 

Figure 2. Spectral energy distributions of the pre-stellar core L1544 and the 
Class protostar IRAS 16293, along with grey-body fits (see text). The 
L1544 SED is based on ISOPHOT, JCMT, and IRAM data from Ward- 
Thompson et al. (1999a,b). The IRAS 16293 SED is based on IRAS, 
ISO-LWS (Correia, Griffin et al. 1999 in prep. - see also Ceccarelli et 
al. 1998), and JCMT data (Sandell 1994). Note that a simple grey-body 
model cannot account for the 25 /jm emission of IRAS 16293 and that a 
two-component model is required (Correia et al. 1999). 

Figure 3. (a) (left) Radial intensity profile of L1689B at 1.3 mm illustrating 
that pre-stellar cores have flat inner density profiles (from Andre, Ward- 
Thompson, & Motte 1996). For comparison, the dotted curve shows a 
spherical isothermal model with p(r) oc r~ 1,2 for r < 4000 AU and p(r) oc 
r~ 2 for r > 4000 AU. The dash-dotted curve shows a model with p(r) oc 
r~ 2 such as a singular isothermal sphere (SIS). 

(b) (right) Typical density profiles expected for a magnetically supported 
core undergoing ambipolar diffusion at various times increasing from to 
to t(j (from Ciolek & Mouschovias 1994). The normalization values are 
n n ,c0 = 2.6 x 10 3 cm' 3 and Rq — 4.29 pc. Open circles mark the instanta- 
neous radius of the uniform-density central region; starred circles mark the 
radius of the magnetically supercritical region (present only for t > fa). 

Figure 4. Plot of lifetime vs. mean density for six core samples compared 
with t oc p -0 ' 75 and t oc p -0 ' 5 as predicted by ambipolar diffusion models 
with different ionisation mechanisms (Jessop & Ward- Thompson 1999). 

Figure 5. Dust continuum mosaic of the p Oph cloud taken at 1.3 mm with 
the IRAM 30m telescope and the MPIfR bolometer array (Motte, Andre, 
& Neri 1998). 

Figure 6. Mass spectrum of the 59 pre-stellar fragments extracted from the 
p Oph 1.3 mm continuum mosaic shown in Fig. 5 (from Motte et al. 1998). 
For comparison, dotted and long-dashed lines show power laws of the form 
AN/AM oc M -1 ' 5 and AN/AM oc M~ 2,5 , respectively. This pre-stellar 
mass spectrum is remarkable in that it resembles the shape of the IMF. 

Figure 7. (a) L bo i-T bo i diagram for 14 well-documented YSOs along with 
three model evolutionary tracks for various (final) stellar masses and cloud 
temperatures (from Myers et al. 1998). Four times t (Myr) since the start 
of infall are indicated, at log t = —1.5, —1.0, —0.5, and 0.0. 
(b) Menv-Lboi diagram for a sample of Class I (filled circles) and con- 
firmed Class sources (open circles) mainly in Ophiuchus, Perseus, and 
Orion (adapted from AM94 and Saraceno et al. 1996a). Evolutionary 
tracks, computed assuming protostars form from bounded condensations 
of finite initial masses and have Ltoi = GMtMacc/R* + £*, where L* is 
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the PMS birthline luminosity (e.g. Stahler 1988), are shown. M env and 
M acc . = M env /r (where r = 10 yr) have been assumed to decline expo- 
nentially with time (see Bontemps et al. 1996a). Small arrows are plotted 
on the tracks every 10 4 yr, big arrows when 50% and 90% of the initial 
condensation has been accreted. The dashed and dotted lines are two M*- 
Lboi relations marking the conceptual border zone between the Class 
{Menv > Mi,) and the Class I (Menv < M*) stage; the dashed line has 
M t oc L bo i (cf. AWB93 and AM94) while the dotted line has M* cx L°^ t 
as suggested by the accretion scenario adopted in the tracks. 
Figure 8. Normalized outflow momentum flux, -Fcoc/Lboi, versus normal- 
ized envelope mass, M env / L°® 1: for a sample of Class (open circles) and 
Class I (filled circles) objects (from Bontemps et al. 1996a). _Fco c/L bo i can 
be taken as an empirical tracer of the accretion rate M acc , while M env /L^ 
is an evolutionary indicator which decreases with time. This diagram 
should therefore mainly reflect the evolution of M acc during the proto- 
stellar phase. The solid curve shows the accretion rate history predicted 
by the simplified collapse model of Henriksen et al. (1997). 
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Table 1: Properties of confirmed Class protostars 



Object 


a(2000) 


5(2000) 


Dist. 








Tb ; 


Outflow 


Infall 


Struc- 


References 








(PC) 


(Lo) 


(M ) 


(%)° 


(K) 


Manifestations 




ture 6 




W30H-TW* 


02:27:04.7 


+61:52:24 


2200 


10 3 -10 4 


~ 20 


~ l 


> 80? 


H2O, radio 


- 


- 


1 


L1448-IRS2 


03:25:22.4 


+30:45:12 


300 


6 


0.9 


3 


70? 


CO 


- 


- 


2 


L1448-N 


03:25:36.3 


+30:45:15 


300 


11 


2.3 


3 


70 


CO, radio 


Y? 


D, B 


3,4,5,6,7,8,9,10 


L1448-C 


03:25:38.8 


+30:44:05 


300 


9 


1.4 


2 


60 


CO, radio, H2 


N? 


D 


3,4,5,11,12,8,10,13 


NGC1333-IRAS2 


03:28:55.4 


+31:14:35 


350 


40 


1.7 


£ 1 


50 


CO, H 2 


Y? 


- 


14,15,9,16,17,13 


SVS13B 


03:29:03.1 


+31:15:52 


350 


~ 7 


2.7 


~ 5 


~ 30 


SiO, H 2 


- 


- 


18,19,20,21,22 


NGC1333-IRAS4A 


03:29:10.3 


+31:13:31 


350 


14 


7. 


5 


34 


CO, radio 


Y 


D, B 


23,24,15,16,25,26,27,28,10,13 


NGC1333-IRAS4B 


03:29:12.0 


+31:13:09 


350 


14 


2.7 


3 


36 


CO, radio 


Y 


D, B 


23,24,15,16,25,28,10,13 


IRAS 03282 


03:31:20.8 


+30:45:31 


300 


1.5 


0.6 


5 


35 


CO, H 2 


- 


- 


29,30,8 


HH211-MM 


03:43:56.8 


+32:00:50 


300 


5 


1.5 


~ 4 


~ 30 


H 2 


- 


- 


31,32 


IRAM 04191 


04:21:56.9 


+ 15:29:46 


140 


0.15 


0.5 


12 


18 


CO, radio 


Y 


N? 


33 


L1527 f 


04:39:53.9 


+26:03:10 


140 


2 


0.4 


0.7 


60 


CO, H 2 


Y 


D 


34,35,36,37,38,6,10,13 


HH114MMS 


05:18:15.2 


+07:12:03 


450 


£ 25 


2.8 


£ 1 


~ 40 


radio, HH 


- 


- 


19 


RN043-MM 


05:32:19.4 


+ 12:49:42 


400 


5 


0.6 


~ 5 


36 


CO, radio, H 2 


N? 


- 


39,19,40,41,10 


OMC3-MM6 


05:35:23.5 


-05:01:32 


450 


< 60 


12 


£ 2 


~ 30 


CO, H 2 


- 


- 


42,43 


L1641-VLAl f 


05:36:22.8 


-06:46:07 


450 


50 


6.5 


~ 3 


70? 


radio, HH 


- 


- 


44,19,45 


NGC2023-MM1 


05:41:24.8 


-02:18:09 


450 


~ 8 


1.8-4.6 


~ 3-10 


~ 30 


CO 


- 


- 


46,47 


NGC2024-FIR5 


05:41:44.5 


-01:55:43 


450 


£ 10 


15. 


£ 30 


20 


CO 


N? 


D, B? 


48,49,50,10 


NGC2024-FIR6 


05:41:45.2 


-01:56:05 


450 


> 15 


6. 


< 10 


20 


CO 


N? 


- 


48,49,50,10 


HH212-MM 


05:43:51.5 


-01:02:52 


400 


14 


1.2 


~ 2 


70? 


CO, H 2 


- 


- 


39,19,51 


HH24MMS 


05:46:08.3 


-00:10:42 


450 


5 


4 


10 


20? 


radio, H 2 


- 


D 


52,53,54,55,9,56 


HH25MMS t 


05:46:07.5 


-00:13:36 


450 


6 


0.5 


5 


34 


CO, radio, H 2 


Y? 


- 


54,57,10 


NGC2264G-VLA2 


06:41:11.1 


+09:55:59 


800 


12 


2 


2 


25 


CO, radio 






58,59,60 


IRAS 08076 f 


08:09:32.8 


-36:05:00 


400 


17 


2.3 


£ 1 


74 


CO, H 2 






61,16,62 


BHR71-MM 


12:01:36.3 


-65:08:44 


200 


10 


2.4 


~ 3 


56 


CO, HH 






63,62 


IRAS 13036 f 


13:07:36.1 


-77:00:05 


200 


1.7 


0.3 


£ 2 


60 


CO 


Y 




62,64,13 


VLA 1623 


16:26:26.4 


-24:24:30 


160 


1 


0.7 


10 


< 35 


CO, radio, H 2 


Y? 


D, B? 


65,66,67,68,69,70,9,10,13,71,22 


IRAS 16293 


16:32:22.7 


-24:28:32 


160 


23 


2.3 


2 


43 


CO, radio 


Y 


D, B 


72,73,74,75,26,27,76,77,10,13 


Trifid-TC3 { 


18:02:07.0 


-23:05:11 


1680 


~ 10 3 


~ 60 


~ 1 


~ 30 


SiO 






78 


L483-MM f 


18:17:29.8 


-04:39:38 


200 


9 


0.3 


~ 0.7 


50 


CO, radio 


Y? 




79,37,10,13 


Serp-S68N 


18:29:47.9 


+01:16:46 


310 


6 


1.0 


3 


40? 


CO, CS 


Y? 




80,81,82,83,84,13 


Serp-FIRSl 1 


18:29:49.9 


+01:15:20 


310 


46 


3 


~ 1 


51 


CO, radio 


N? 


N? 


85,80,81,82,86,84,10,13,87 



Table 1 (cont'd) 



Object 


a(2000) 


5(2000) 


Dist. 






^-'smm l^bol 


Tb ; 


Outflow 


Infall 


Struc- 


References 








(PC) 


(Lo) 


(M ) 


(%)° 


(K) 


Manifestations 




ture 6 




Serp-SMM4 


18:29:57.1 


+01:13:15 


310 


9 


3 


~ 3 


35 


CO, radio 


Y 


N? 


85,81,82,84,88,10,13,87 


Serp-SMM3 f 


18:29:59.7 


+01:14:00 


310 


8 


0.9 


~ 1 


40 


CO, H 2 


N? 


N? 


85,81,82,84,89,10,87 


G34.24+0.13MM* 


18:53:21.5 


+01:13:45 


3700 


4000 


100 


~ 0.5 


50 


CH 3 OH 






90 


L723-MM 


19:17:53.7 


+ 19:12:20 


300 


3 


0.6 


~ 4 


50 


CO, radio 


N? 




91,92,32,19,93,94,10 


B335 


19:37:00.8 


+07:34:11 


250 


3 


0.8 


6 


37 


CO, radio 


Y 


D? 


91,95,96,40,97,10,13,98 


S106-SMM 


20:27:25.3 


+37:22:46 


600 


£ 24 


£ 10 


£8 


£ 20 


bip. HII, H 2 






99,70 


L1157-MM 


20:39:06.2 


+68:02:22 


440 


11 


0.5 


~ 5 


~ 60 


CO, H 2 


Y 


D? 


100,32,101,69,10,13 


GF9-2 


20:51:30.1 


+60:18:39 


200 


0.3 


~ 0.5 


~ 10 


£ 20 




Y 




102,103,104 


CepE-MM 


23:03:13.1 


+61:42:26 


730 


75 


7 


~ 2 


~ 60 


CO, H 2 






105,106 


IRAS 23385* 


23:40:54.5 


+61:10:28 


4900 


16000 


370 


~ 0.7 


£ 40 


SiO 






107 



L smm is the luminosity radiated longward of 350 (im; Class 0s are defined by L smm /L(, ; > 0.5% (see III— A). 

6 Small-scale structure: D = presence of disk-like component; B = binary or multiple system; N = single object with no disk 

t Border-line Class 0. * Candidate massive Class object. 
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